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Abstract 

163 Ho to 163 Dy is probably due to the small Q value of about 2.5 keV the best case to determine 
the neutrino mass by electron capture. The energy of the Q value is distributed between the 
excitation of Dysprosium (and the neglected small recoil of Holmium) and the relativistic energy 
of the emitted neutrino including the rest mass. The reduction of the upper end of the deexcitation 
spectrum of Dysprosium below the Q value allows to determine the neutrino mass. The excitation 
of Dysprosium can be calculated in the sudden approximation of the overlap of the electron wave 
functions of Holmium minus the captured electron and one- , two-, three- and multiple hole- 
excitations in Dysprosium. Robertson [R. G. H. Robertson, Phys. Rev. C91, 035504 (2015) and 
arViv: 1411.2906] and Faessler and Simkovic [Amand Faessler, Fedor Simkovic, accepted for Phys. 
Rev. C, March 2015 and arXiv: 1501.04338] have calculated the influence of the two-hole states 
on the Dysprosium spectrum. Here for the first time the influence of the three-hole states on the 
deexcitation bolometer spectrum of 163 Dysprosium is presented. The electron wave functions 
and the overlaps are calculated selfconsistently in a fully relativistic and antisymmetrized Dirac- 
Hartree-Fock approach in Holmium and in Dysprosium. The electron orbitals in Dy are determined 
including the one-hole states in the selfconsistent iteration. The influence of the three-hole states 
on the Dy deexcitation (by X-rays and Auger electrons) can hardly be seen. The three-hole states 
seem not to be relevant for the determination of the electron neutrino mass. 


1 Introduction 

One of the most pressing problems in particle physics is the determination of the absolute value of the 
neutrino mass. The single beta decay, specifically the Triton decay, can give the electron antineutrino 
mass P (present upper limit [2]: rn 9 < 2.2 eV), while the neutrinoless double beta decay yields 
the mass of the effective Majorana electron neutrino mass [3], if the transition matrix element can be 


calculated reliably. At the moment the upper limit for the Majorana neutrino mass from the double 
beta is about 0.3 eV. 

The electron neutrino mass can be determined by electron capture. The most favorable case due 
to the small Q value of about 2.5 keV seems to be capture in 163 Ho (Z=67) to atomic excited states 
in 163 Dy (Z=66). The energy of this Q value is divided into the excitation of the Dysprosium atom, 
(plus the recoil of the Holmium atom, which can be neglected,) and into the energy carried away by 
the electron neutrino (kinetic energy and rest mass). The deexcitation spectra (X rays and Auger 
electrons) of all excited states in the Dysprosium atom have the upper end at the Q value minus 
the neutrino rest mass. The determination of the electron neutrino mass from electron capture is 
more complicated than the determination of the antineutrino mass from the Triton decay, since more 
parameters have to be fitted to the data at the upper end of the spectrum. Under the assumption, that 
one excited Dy resonance with a Lorentzian line shape determines the behavior at the upper end, one 
has to fit simultaneously four parameters to the experiment after the detector response is folded into 
the theoretical spectrum: (1) The neutrino mass, (2) the energy distance of the leading resonance to 
the Q value, (3) the width and (4) the strength of the resonance. Originally only one-hole excitations 
in the Dysprosium atom have been included in the theoretical determination of the spectrum HIS- 
Recently Robertson [6] and with major improvements Faessler and Simkovic [7j have included two-hole 
excitations in Dysprosium. The main improvement of ref. [7j over ref. [6] was, that the electron wave 
functions for Holmium and for Dysprosium have been calculated in these nuclei (Z=67 and Z=66) 
selfconsistently in a fully relativistic and antisymmetrized Dirac-Hartree-Fock approach mum and 
not in Xenon (Z=54) as by Carlson and Nestor [11], [12j , results used by Robertson [6|. Due to the 
additional occupied states in Holmium and Dysprosium more two-hole states are allowed. In addition 
the electron orbitals in Dysprosium are determined selfconsistently with the appropriate hole [7]. The 
finite size of the nuclear charge distribution is included using the Fermi parametrization determined by 
electron-nucleus scattering. The theoretical formulation was derived in second quantization including 
automatically the antisynmretrization and the overlap and exchange corrections. 

In the present work this formulation is extended for the first time to the three-hole states in Dys¬ 
prosium including again all the improvements of the description of the two-hole states of ref. [7]. 


2 Description of Electron Capture and 
the Excitation of the Three-Hole States. 

The calorimetric spectrum of the deexcitation of 163 Dy after electron capture in 163 Ho can be expressed 
according to refs. |T] and m-- 
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The second sum over the excited states h’ in Dy includes the one-hole states b’, the two-hole 
states b'{j)~ l q ) and the three-hole states b'(p^, q v )(p^ 1 qy)■ With the Q value Q = 2.3 to 2.8 keV 





unmans eei mi and a recommended value pH] Q = (2.55 ±0.016) keV. The ECHo collaboration [US] 
measured for the Q value for electron capture in 163 Ho to ie3 Dy: 

Q(ECHo) = 2.80 ± 0.08 keV. (3) 

t/g j is the probability for the admixture of different neutrino mass eigenstates i = 1, ,.N V into the 
electron neutrino and E c is the excitation energy of final Dysprosium. B b with the same orbital quantum 
numbers b = b’ in Ho and Dy are the overlap and exchange corrections [5]; Ao contains the nuclear 
matrix element squared [21]; E h > are the one-, two- and three-hole excitation energies in Dysprosium. 
T h > is their widths [5]. 

We assume, that the total atomic wave function can be described by a single Slater determinant. B b 
takes into account the overlap and the exchange terms between the parent \G > and the daughter atom 
in the state | A' b , > with a hole in the electron state | b' >. We use the sudden approximation as Faessler 
et al. [22] . B b for the electron capture probability from the state b relative to the capture from 3 si/ 2 
with one hole in b’ in the Dysprosium atom is given in eq. (jTJ) in the Vatai approximation [23, 21] ■ But 
the numerical value used here are calculated with the full overlap and exchange corrections of Faessler 
et al. [5]. 
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For two-hole final states one has to multiply eq. (|4]) with the probability to form a second hole 
characterized by the quantum numbers ”p”. One has to replace < A' b ,\ai\G > by 

< A' b , , ,\cii\G > with the two electron holes b’ and p’ and the additional electron particle q’ in 
Dysprosium above the Fermi surface F. These expressions are derived in reference |7j. The formulas for 
three-hole states can be obtained with the help of Wick’s theorem [25j from: 
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With the Holmium ground state: 


(5) 


|G >= a{a2c4...G4|0 > (6) 

The probability for forming three hole states can be calculated in the sudden approximation with 
the help of Wick’s theorem [[25] . 

To evaluate the probability for two electron particle-hole states one sums incoherently over all 
unoccupied states q\ and q’ 2 ■ 
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Here as stressed above the sum over q[ and q' 2 runs over the unoccupied bound and continuum states 
in Dy. The dependence on the quantum numbers b = b’ originates from the fact, that the electron 
orbitals in Dysprosium are calculated selfconsistently with an empty state b’. One can now use the 




completeness relation to shift the sum over q[ from the continuum to states, which one can calculate 
easier. One divides the completeness relation into two pieces: up to the last occupied state below the 
Fermi Surface F and all states above the last occupied state including also the continuum. 
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The sum in eq. (J7J) is the last part of the completeness relation (JHJ) and one can transcribe (J7J) into: 
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An interchange of two electrons from the states b', p\ and p' 2 yields due to the antisymmetrization 
a sign. But for the probability the minus sign is irrelevant. In the literature one uses often the 
Vatai approximation [231 [24|: Exchange corrections have been neglected already in eq. (]9|) . In addition 
one assumes, that the overlaps of electron wave functions in the parent and the daughter atom with 
the same quantum numbers can be approximated by unity. Typically the overlaps have values [5] 
< k ’ Dy \ k Ho >~ 0.999 and 0.999 65 ~ 0.94. In the Vatai approximation |24] one replaces this value 
by 1.0. The probability for two one-electron particle and one-electron hole states is the product. The 
probability for a second hole in p[ with a first hole in b’ is: 
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The physics of the two terms subtracted from 1 in eq. (fTO is: The Erst subtracted term gives the 
probability, that the state p\ in Dy is occupied. The sum in the second term takes into account the Pauli 
principle and prevents, that electrons can be moved into occupied states in Dy. The single electron 
states like |p >= \n,£,j,m > include also the angular momentum projection quantum number m. This 
projection is for the description of the data irrelevant. The Erst subtracted term in eq. (fTO gives the 
probability, that a specific magnetic substate m[ is occupied in p t \ The probability, that all magnetic 
substates of pi are occupied, is the product of the probabilities of all substates. One obtains the Nth. 
power of the single electron probability with N p / = N n ^j. p ', = (2 j + 1 ) p r. 
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for n and n’ with: N n j^ = N n i^j = 2 j + 1; 


With the expression: 


(12) 



For the hole state in b’ and the partially occupied state in 4 f 7 / 2 < one needs special expressions: 

N( n ,e,j) b = 2 Jb and N 4f7 / 2 < = 5. (13) 

Here \b > is the orbital in Ho, from which the electron is originally captured. The electron p\ is moved 
to 4/7/2 with now 5 electrons in this orbit. Thus the second (and the third) hole is in Dy in the p\ 
electron orbit. N n ^j /(2 j +1) is the averaged probability to find an electron in the | p, n, £, j, m > orbital 
and N n r^j the number of electron in the \q' i ,n',£,j > state. Primes always indicate states in Dy. For 
the probability of the second and third hole we used the so called Vatai |23f [2j] approximation with the 
overlaps of corresponding electron wave function with the same quantum numbers in Ho and Dy equal 
to unity and neglected the exchange corrections. For the ’’diagonal” overlaps of the order of ~ 0.999 
this is a good approximation. 

For the probability of two different electron-particle electron-hole excitations <Zi(pi) _1 and q 2 (p 2 )~ 1 
in eq. (J7]), eq. ([9]) one has now to introduce eq. (TTOl) or finally eq. (TTIT) as product for particle 1 and 
particle 2 into eq. ((9]). 
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The three hole probability normalized to the electron capture from the Ml 3sl/2 state in 163 Holmium 
is then obtained by multiplying this expression into eq. (HJ) for P&. 


3 Spectrum of Dy deexcitation with three-hole states. 

Tables [Tj and [2] show the one and the two-hole resonance energies and widths taken from results of the 
ECHo collaboration [161 T9] and, if there not available, taken in reference [6] from the handbook of 
Chemistry tabulated by Weast [26]. The two hole energies are determined as the one hole excitation 
energy in Dysprosium and the electron-particle and electron-hole energy differences in Holmium. Due 
to the hole in an inner shell in Dysprosium, which reduces the shielding of the charge of the nucleus, 
the field looks for the outer electron particle-hole excitation in Dysprosium similar as in Holmium with 
one positive charge more. This is only a rough approximation, since the relevant one-hole states in 
Dysprosium Ml 3sl/2, N2 3pl/2, N1 4sl/2, N2 4pl/2, 01 5sl/2 and 02 5pl/2 are not from the 
inner most shells and thus the outer particle-hole states in Dysprosium do feel only approximately the 
selfconsistent field of Holmium. A better approach would be to calculate the electron particle-hole 
energies in Dysprosium selfconsistently with the hole in the appropriate state. This could shift the 
resonance energies by a few eV compared to the approximation used. 

The three hole resonance energy is in the approximation discussed: 




Table 1: The 6 one-hole excitation energies Ec and the widths T in 163 Dy with quantum numbers 
n, £, j are listed according to the ECHo collaboration [13] [16] [191HQ] as far as they are available. Else 
they are taken according to ref. [6j from Weast tabulated in the Handbook of Chemistry [26], although 
there seem to be in some cases better values in the literature [29], [30], |3Tj and [32]. A Q value of 
2.8 keV allows energetically the excitations of the one hole states: Ml 3sl/2, M2 3pl/2, N1 4sl/2, N2 
4pl/2, 01 5sl/2 and 02 5pl/2. Electrons can be captured only from sl/2 and pl/2 orbits in Holmium, 
since other electrons have no probability to be at the nucleus. 


Number 

1. hole 

E c [eV} 

r[eV] 

Number 

1. hole 

E c 

r[eV] 

1 

3sl/2 

2040 

13.7 

2 

3pl/2 

1836 

7.2 

3 

4sl/2 

411 

5.3 

4 

Apl/2 

333 

8 

5 

5sl/2 

48 

4.3 

6 

hpl/2 

30.8 
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Where the particle-hole energy can be taken from table [2] starting with column four and subtracting 
the corresponding one hole energy E cb -i from table [1] 

E Ct particle—hole energy \Ec,(pi)~ 1 ,Ho ^c,q\,Ho) 

Ecjb.jn)- 1 (column 4, table \ 2]) — E cb -i(one — hole energy of tabled (16) 

The numerical calculation includes 6 Dysprosium one-hole states 3s 1/2, 3pl/2, 4sl/2, 4pl/2, 5sl/2 
and bpl/2 listed in table [Tj and 39 two-hole states tabulated in table [2] The problem of choosing for 
the particle-hole energies the values of Ho is discussed before eq. (1131) . The allowed 208 three hole 
states in this work are limited by energy conservation to Q > Eq( 1 — hole) + Ec(particle — hole i) + 
Eq (particle—hole 2 ) with the particle-hole excitation energies Ec(particle—hole ), taken as the ionization 
energy in Ho. 


Figure [U shows the sum of the one-, two- and three-hole theoretical bolometer spectrum after electron 
capture in Holmium for the deexcitation of Dysprosium. The electron wave functions are calculated 
in a Dirac-Hartree-Fock approach la la mg. In Dysprosium the one-hole states are included in the 
selfconsistent determination of the electron wave functions. A detailed description of our approach is 
given in ref. [7|. Figure [2] compares the sum of the theoretical one- plus two- plus three-hole deexcitation 
spectrum of Dysprosium with the results including only one-hole states and the sum of the one- and 
the two-hole states. The three-hole contributions can practically not be seen even in this logarithmic 10 
plot. 

Figure [3] compares in a logarithmic 10 plot the deexcitation spectrum of the one-hole states, the 
one- plus two-hole and the one- plus two- plus three-hole excitations with the ECHo data [191, 20] . The 
experimental spectrum in binned into 2 eV. Since some bins contain no counts the logarithmic value is 
— 00 . For a more relevant logarithmic plot one would need more counts. The three experimental peaks 
between 1.2 and 1.6 keV originate from a Ui Pm contamination. Figure [4] shows in detail a comparison 
of the one-, the two- and three-hole spectrum with the ECHo data HSIEQ] around the area of the N2 
4pl/2, 0.333 keV and the N1 4sl/2, 411 keV hole states. 
















Table 2: The 39 (from number 7 to 45) two-hole excitation energies Eq and the widths T in 163 Dy 
with quantum numbers n, £, j are listed. We adopt here the first hole values for E c and the width T 
of the ECHo collaboration [131 EES IS [2D] and add the excitation energies of the second hole taken by 
Robertson [6] taken from Weast tabulated in the Handbook of Chemistry [26]. The data measured by 
the ECHo collaboration [19] , [20] are listed in table [0 The energies of the two hole states are given by 
Ec(lholeDy ) + Ec(particle — hole excitation in Ho). The number of one- and two-hole states listed 
in tables |T] and [2] is 45. The allowed 208 three hole states included in this work are limited by energy 
conservation to Q > E c (l — hole ) + E c (particle — hole i) + E c (particle — hole 2 ). The particle-hole 
excitation energies Ec(particle — hole ) can be obtained from this table by subtracting the one-hole 
excitations given in table [0 


Number 

1. hole 

2. hole 

E c [eV] 

r[eV] 

Number 

1. hole 

2. hole 

E c 

r[eV] 

7 

3sl/2 

4sl/2 

2472.4 

13.7 

27 

4sl/2 

4sl/2 

841.4 

5.4 

8 

3sl/2 

4pl/2 

2385.3 

13.2 

28 

4sl/2 

4pl/2 

752.5 

5.4 

9 

3sl/2 

4p3/2 

2350.0 

13.2 

29 

4sl/2 

4p3/2 

717.2 

5.4 

10 

3sl/2 

4d3/2 

2201.8 

13.2 

30 

4sl/2 

4d3/2 

569.0 

5.4 

11 

3sl/2 

4d5/2 

2201.8 

13.2 

31 

4sl/2 

4d5/2 

569.0 

5.4 

12 

3sl/2 

4f5/2 

2050.4 

13.2 

32 

4sl/2 

4f5/2 

417.6 

5.4 

13 

3sl/2 

4f7/2 

2047.0 

13.2 

33 

4sl/2 

4f7/2 

414.2 

5.4 

14 

3sl/2 

5sl/2 

2091.1 

13.2 

34 

4sl/2 

5sl/2 

458.3 

5.4 

15 

3sl/2 

5pl/2 

2072.6 

13.2 

35 

4sl/2 

5pl/2 

439.8 

5.4 

16 

3sl/2 

5p3/2 

2065.9 

13.2 

36 

4sl/2 

5p3/2 

433.1 

5.4 

17 

3pl/2 

4sl/2 

2269.2 

6 

37 

4pl/2 

4pl/2 

671.8 

5.3 

18 

3pl/2 

4pl/2 

2180.3 

6 

38 

4pl/2 

4p3/2 

636.5 

5.3 

19 

3pl/2 

4p3/2 

2145.0 

6 

39 

4pl/2 

4d3/2 

488.3 

5.3 

20 

3pl/2 

4d3/2 

1996.8 

6 

40 

4pl/2 

4d5/2 

488.3 

5.3 

21 

3pl/2 

4d5/2 

1996.8 

6 

41 

4pl/2 

4f5/2 

336.9 

5.3 

22 

3pl/2 

4f5/2 

1845.4 

6 

42 

4pl/2 

4f7/2 

328.3 

5.3 

23 

3pl/2 

4f7/2 

1842.0 

6 

43 

4pl/2 

5sl/2 

377.6 

5.3 

24 

3pl/2 

5sl/2 

1886.1 

6 

44 

4pl/2 

5pl/2 

359.1 

5.3 

25 

3pl/2 

5pl/2 

1887.6 

6 

45 

pl/2 

5p3/2 

352.4 

5.3 

26 

3pl/2 

5p3/2 

1860.9 

6 

— 

— 
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Figure 1: (Color online) Logarithmic 10 three-hole Bolometer Spectrum (J2]) with (J4| and 
(fT4l) for the sum of the one-, the two- and the three-hole probabilities calculated in this work 
with the resonance energies and widths of tables Q] and [2] with the assumed Q-value Q = 
2.8 keV for the bolometer energy between 0.0 and 2.8 keV. The logarithmicio coordinates 
of the ordinate has to be read as io ordmate , e. g. the ordinate value -2 corresponds to 10 -2 . 
The theoretical spectra are normalized to the experimental 4 si/ 2 hole peak at 0.411 keV (see 
figure ED). 
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Figure 2: (Color online) Logarithmic^ plot of the bolometer spectrum (J2J) including the one- 
, the one- plus two- and the one- plus two- plus three-hole states. The Q value is assumed 
to be Q = 2.8 keV according to an ECHo measurement [TB]. The resonance energies and 
widths are listed in tables [T] and [2j The three-hole states can hardly be seen. 
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Figure 3: (Color online) Comparison of the logarithmiciQ spectrum of the one-, the one- plus 
two and the one- plus two- plus three-hole states with the data of the ECHo collaboration 
mmmm- The experimental data are binned in 2 eV (counts/2eV). Some bins have no 
counts and thus the logarithmic value is — oo. The effect of the three-hole states does show 
hardly up even in this logarithmic figure. 
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Figure 4: (Color online) Comparison of the logarithmicio bolometer spectrum around the 
N 1 = (4sl/2) -1 0.411 keV (relative probability to the 3sl/2 state 24.4 %) and N2 = 
(4pl/2) -1 0.333 keV (relative probability 1.22 % ) resonances with the ECHo data [13l lT6lH9l 
I20jj. At 0.569 keV one sees the degenerate two-hole states IV1 = (dsl/2)^ 1 ,1V4 = (4d3/2) -1 
with the relative probability of 0.088 % and N 1 = (Asl/2)~ l , N5 = (4d5/2)^ 1 with the 
relative probability of 0.125 % (total 0.213 %). The resonance energies and the widths are 
taken from the ECHo collaboration [T9l '20j • The theory does not include the hnite resolution 
of the measurement, which is about 10 eV FWHM. The sum of the one- plus two- and the 
sum of the one- plus two- plus three-hole spectra look even in a logarithmic scale identical. 











4 Conclusions 


Here for the first time the effect of three hole excitations on the deexcitation spectrum of 163 Dy after 
electron capture in 163 Ho has been studied. The spectrum is not affected visibly by the three-hole 
excitations of electron hole and two electron particle-hole states. The two electron particle-hole exci¬ 
tations plus the energy of a one-hole state (three-hole states) add up to higher energies as the one- 
and two-hole states. A one-, two- or three-hole state close to the Q value could complicate the deter¬ 
mination of the neutrino mass. If only one resonance state is dominant near the Q value and if the 
line profile is Lorentzian (or of an other known analytical form with two parameters), one needs to fit 
simultaneously four parameters of the theory to the data: (1) neutrino mass, (2) the energy difference 
of the dominant state to the Q value, (3) the width of the resonance and (4) the strength. Before fitting 
the four parameters of the theoretical spectrum to the data, one has to fold the experimental spectral 
function of the detector into the theoretical results. The situation is more complicated, if several one- 
or two- or three-hole resonances contribute to the shape of the spectrum near the Q value. With the 
recent measurement [19, 20] of the ECHo collaboration (J3]) Q = 2.80 ± 0.08 [fceK], it has been shown 
by Faessler and Simkovic [7], that the highest one hole state is dominant at the Q value. The high 
energy two-hole states do not influence the behavior in the region, which determines the neutrino mass 

[7]. Although the three-hole excitations can come closer to the Q value, their strength is so weak, that 
they are not dangerous for the determination of the neutrino mass. 


Acknowledgment: We want to thank other members of the ECHo collaboration for discussions. 
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